Instabilities in the Bogoliubov Spectrum of a condensate in a 1-D periodic potential 
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We study the stability of standing wave solutions to a one- 
dimensional Gross-Pitaevsky equation with a periodic poten- 
tial. We use some simple complex analysis and the Hamilto- 
nian structure of the problem to give a simple rigorous cri- 
terion which guarantees the existence of non-real spectrum, 
which corresponds to exponential instability of the standing 
wave solution. This criterion can be stated simply in terms 
of the spectrum of one of these self-adjoint operators. When 
the standing wave has small amplitude this criterion simpli- 
fies further, and agrees with arguments based on the effective 
mass in the periodic potential. 



The possibility of creating cold atom condensates in 
optical lattices has made it possible to observe both linear 
phenomena, such as Bloch oscillations [1] as well as gen- 
uinely nonlinear phenomena such as soliton generation 
[2] and an insulating transition [3]. While a great deal 
of effort has gone into the construction of exact or ap- 
proximate solutions the question of stability is foremost, 
since it dictates what states can actually be observed, 
as well as governing the evolution of small disturbances 
to these states. Recent papers of Wu and Niu [4,5] and 
Burger et.al. [6] have helped elucidate the role of the dy- 
namic and Landau instabilities. These dynamic instabil- 
ities have also been realized experimentally [7]. Most of 
the current theoretical stability work is either numerical 
or relies on either special properties of some exact solu- 
tion [8-10], or some approximation, most notably a two- 
level approximation [5], the tight binding approximation 
[5,11,12], weak nonlinearity [13], or long wavelength [14]. 
In this paper we give a simple, rigorous, easily checked 
sufficient condition for the existence of exponential insta- 
bilities. In the case of a weak nonlinearity this condition 
reduces to a physically natural condition in terms of the 
effective mass of the lattice. 

Our starting point is the 1-D Gross-Pitaevsky equation 
with a periodic potential, 



-~d xx + V(x) + \(f>\ 2 (x) -u))v = L+v 



iipt = 



±\i}\ 2 ip + V(x)ip V(x + 1) = V(x) 



with a standing wave solution of the form ip(x,t) — 
exp(—iojt)(f>(x), where the standing wave profile <j> is real 
and either periodic 4>{x + 1) = <fc(x) or antiperiodic 
4>(x + 1) = —cj)(x). The linearized equation governing 
small disturbances is 



iv t = I- -d xx + V(x) + 3\(t)\ 2 (x) - luJu = L_u 

where u, v are the real and imaginary parts of the com- 
plex disturbance. Upon taking a Fourier transform in 
time this takes the form of the well-known Bogoliubov 
equation 



[XU 
(XV 



Tj + V 

L_u. 



Because of the Hamiltonian structure of the problem, the 
above eigenvalue problem can be written in vector form 

v x = JB.(x, li)v 

with v = (u,v,u x ,v x Y , the matrix H = is Hermi- 
tian, and the matrix J is the canonical Hamiltonian form 
I 
-I 

guarantee that fx is in the spectrum of the operator above 
if and only if the period map M(/i) has eigenvalues on 
the unit circle, where the period map M(/x) is defined by 



The standard results of Floquet theory 



U x =JH(x, M )U U(0, M ) = I 
M(//) = U(1, M ) 



(1) 
(2) 



In the second order case, which arises in the classical 
Floquct-Bloch theory, the location of the eigenvalues of 
the period map are determined entirely by the trace of the 
(2 x 2) matrix M, there being two eigenvalues on the unit 
circle if Tr(M) £ (—2, 2) and two eigenvalues off of the 
unit circle if Tr(M) ^ [-2,2]. The case of the Bogoliubov 
operator is somewhat more complicated, since it consists 
of coupled second order problems, but can be treated 
in much the same way due to the Hamiltonian structure 
[15]. First we note that U and thus M satisfy the relation 

U t (x,Li)3V(x,Li) = J, 

a result known as the Poincare-Liouville theorem. From 
this it follows that M is similar to its inverse transpose. 
This implies that the eigenvalues are invariant under in- 
version with respect to the unit circle, which in turn im- 
plies that the characteristic polynomial of M has the fol- 
lowing form: 



det(M - AI) 
a 



l + a\ + b\ 2 + a\ 3 + A 4 

- Tr(M) 

- (Tr(M) 2 - Tr(M 2 )) . 



(3) 
(4) 

(5) 
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Because of the special form the above, the quartic admits 
an explicit factorization into quadratics, 

det(M - AI) = (1 - K+(fi)X + A 2 )(l - A + A 2 ), 

where the quantities K± are given by 

Tr(M) ± i/2TrM 2 (/i) - ( Tr M^)) 2 + 8 
K±{p) = g • 

It then follows that the period map for the Bogoliubov 
operator has two eigenvalues on the unit circle if K + is 
real and K + e (—2, 2), and another two on the unit circle 
if K_ is real and K_ e (—2,2). It can be shown using 
complex analysis that the spectrum of the Bogoliubov 
operator consists of a union of piecewise smooth arcs in 
the complex plane, where the endpoints of the arcs are ei- 
ther band edges K±(fi) — ±2, critical points K'±(fj,) = 0, 
or branch points K + (ii) = K-(p). It is worth noting 
that a generic quartic does not admit a factorization into 
quadratics of the above form. It is only because of the 
Hamiltonian nature of the problem that this occurs. For 
a fourth order problem without this symmetry it is nec- 
essary to use the general solution to the quartic, and a 
calculation analogous to the present would be exceedingly 
difficult. 

We are primarily interested in instabilities of modu- 
lational type, where the Bogoliubov operator has spec- 
trum off of the real axis in a neighborhood of \i = 0. It 
is also possible for the operator to exhibit sideband type 
instabilities, but these instabilities are more difficult to 
detect. We will consider these in a future paper. To de- 
rive a criterion for a modulational instability, we note the 
following facts 

• Tr(M 2 (^i)) , Tr(M(^)) 2 are even, analytic func- 
tions of a complex fi 

• K±{p) are even and analytic in a neighborhood of 
A* = if K+(0) ^K-(0). 

• L_ (considered as a periodic Schrodinger operator) 
has a band edge at /i = 0. 

The first follows from the fact that the Bogoliubov op- 
erator is invariant under the transformation u — > —u, v — > 
v,fi — > — (i. This implies that the monodromy ma- 
trix M(/z) satisfies M(/i) = VM(-/i)V* with V or- 
thogonal, which implies that Tr(M(/i)) = Tr(M(- / u)) 
and likewise for Tr ( M 2 (/i)) . This symmetry reflects 
the symmetry between left and right propagating waves. 
This symmetry is not present if one considers the lin- 
earization about topological-type standing waves which 
possess a non-trivial quasi-momentum, and the analy- 
sis is much more difficult. Standard arguments from 
the theory of ordinary differential equations show that 
Tr(M(/x)) is an analytic (matrix-valued) function of a 
complex \i. It follows from the formula for K± that, 



away from the branch points where (K + (/i) — K-(n)) 2 — 
2TrM 2 (^) - (TrM(p)) 2 + 8 = 0, the Floquet discrimi- 
nants K±(n) are analytic functions. The final fact follows 
from the U(\) symmetry of the Gross-Pitaevsky equation 
ip — > ipe l9 , which (via Noether's theorem) implies that 
L_4> = 0- Since <f> is either periodic or antiperiodic this 
implies that zero energy is a band edge for L . 
The main result of this paper is the following: 

Theorem 1 Suppose that $ is a standing wave solution 
to the GP equation for which the L + operator is in the 
interior of a band. Then this solution to the GP equation 
is exponentially unstable. 

Proof: When fi = the equations for u and v decou- 
ple, and the period map for the GP equation is block 
diagonal, with one 2x2 block (denoted by niL + corre- 
sponding to L_|_, and the other block (ixil_) correspond- 
ing to L . If fi = is not a band edge for L + it follows 

that (K+(0) - iMO)) 2 - (Tr(m L J - Tr(m L+ )) 2 + 0, 
so the Floquet discriminants are analytic in a neighbor- 
hood of [i = 0. Since K± are even functions it follows 
that both are real on some interval of the imaginary axis. 
If fi = is in the interior of a band for L + it follows that 
.K+(0) G (—2, 2), and so there is some interval along the 
imaginary axis on which K + is real and e (—2, 2). Thus 
the GP equation has a band of spectrum along the imag- 
inary axis. 

An important special case of this is given by the fol- 
lowing: 

Corollary 1 Small amplitude solutions of the focusing 
GP equation corresponding to lower band edges are mod- 
ulationally unstable, small amplitude solutions of the de- 
focusing GP equation corresponding to upper band edges 
are modulationally unstable. 

Proof: The proof is perturbative, considering L + as a 
perturbation of L_. Since L + =L_±2|0| 2 the pertur- 
bation is strictly positive (resp. negative) and it follows 
that the band edges of the L + and L_ operators satisfy 
At„(L+) > /i„(L_) (resp. n„(l>+) < /u n (L_)). Since the 
unperturbed operator L_ has a band edge at fi = it 
follows that, for \\(fi\\oc « 1, M = will be in the inte- 
rior of a band of L + for the focusing case and a lower 
band edge, or defocusing case and an upper band edge. 
From the above result this guarantees the existence of an 
instability. 

This result provides a rigorous justification of the in- 
tuition provided by the effective mass theories, such as 
those proposed by Konotop and Salerno [13] and Tay- 
lor and Zaremba [14] where the effective mass is given 
by the curvature of the dispersion relation. Since this 
curvature has the same sign as the bare mass at lower 
band edges and the opposite sign at upper band edges 
this kind of argument leads to the same conclusion. An- 
other way to look at this result is that in either of the 
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above cases (lower band edge with a focusing nonlinear- 
ity or upper band edge with a defocusing nonlincarity) 
the periodic solution is unstable to formation of a gap 
soliton. The obvious advantage of the more general cri- 
terion is that it applies for strong nonlinearity, when the 
band-gap structure of L± can differ considerably from 
the linear problem, and effective mass arguments based 
on the linear problem can no longer be expected to hold. 

We have conducted some numerical experiments to il- 
lustrate these results. For our basic model we take the 
one dimensional Gross-Pitaevsky equation with a Jacobi 
elliptic function potential: 

iipt = -\^ xx ± HV + V sn 2 (x,k)ip. 

As is shown in previous work [8] this equation has a fam- 
ily of elliptic function solutions. In the first experiment 
we take the focusing sign of the nonlincarity. In this case 
we have an exact solution 



{Re K + , Re K_ } 



{ImK + , ImK_ } 



ip(x, t) = yJ-(V + k 2 )sn(x, k) exp(-iwi) V < -k 2 

which represents a solution bifurcating from the third 
band edge - the lower band edge of the second band. 
The corollary implies that for | Vb + k 2 \ « 1 there ex- 
ists a band of spectrum along the imaginary axis. Fig- 
ure 1 shows a plot of the Floquet discriminants K±(p) 
along the real \x axis for k 2 = .5, Vq = —.6. Note that at 
K_(0) = —2, via Noether's theorem and the fact that the 
third band edge represents an antiperiodic cigenfunction. 
The other Floquet discriminant satisfies ^+(0) w —1.4, 
and is thus in the interior of a band. From the first theo- 
rem it follows that there is a band of spectrum along the 
imaginary \x axis. Figure 2 shows a plot of K±(fi) along 
the imaginary axis. Both Floquet discriminants are real 
and 6 (—2, 2) near the origin, implying that the Bogoli- 
ubov eigenvalue problem has a band of spectrum along 
the imaginary axis. This corresponds to an exponential 
instability in the evolution of small disturbances to this 
standing wave solution. 
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FIG. 1. Plot of the Floquet discriminants K±{li) for /i 
along the real axis for an sn(a;, k) type solution to the focus- 
ing GP equation with an elliptic function potential. Note 
that K+(Q) e (-2,2) 
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FIG. 2. Plot of the Floquet discriminants K±(/j,) for fi 
along the imaginary axis for an sn(x, k) type solution to the 
focusing GP equation with an elliptic function potential. 
Note that K±(li) are real in an interval of the imaginary jx 
axis containing the origin. The arc of spectrum terminates 
at a pair of branch points at [x » ±.07i 

In summary, by using the Hamiltonian structure and 
symmetries of the Gross-Pitaevsky equation with a pe- 
riodic potential, we have presented a sufficient condition 
for the exponential instabilities of standing wave solu- 
tions. 

It is our hope, that this criterion will be a complement 
to the numerous already known solutions in studying the 
Gross-Pitaevsky equation, and condensates in optical lat- 
tices. 
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